Glassy Ratchets For Collectively Interacting Particles 



cn 

O 

o 
< 

in 



O 



I 

o 



ON 

o 

00 

o 
m 
o 



o 



C. Reichhardt, C.J. Olson Rcichhardt, and M.B. Hastings 
Center for Nonlinear Studies and Theoretical Physics Division, Los Alamos National Laboratory, Los Alamos, NM 87545 

(February 2, 2008) 

We show that ratchet efTects can occur in a glassy media of interacting particles where there is no 
quenched substrate. We consider simulations of a disordered binary assembly of colloids in which 
only one species responds to a drive. We apply an asymmetric ac drive that would produce no net 
dc drift of an isolated overdamped particle. When interacting particles are present, the asymmetric 
ac drive produces a ratchet effect. A simple model captures many of the results from simulations, 
including flux reversals as a function of density and temperature. 
PACS numbers; 82.70.Dd, 05.40.-a, 05.70.Ln 



In ratcheting systems, a dc current or flux arises upon 
application of a periodic drive, such as when the underly- 
ing potential substrate is flashed, or, for a flxed substrate, 
when alternating forces are applied [1]. Ratchet effects, 
which can be deterministic or stochastic in nature, have 
been studied in a wide variety of systems, including bi- 
ological motors [2], flux motion in superconductors [3], 
granular transport [4], and colloidal transport. In many 
of these systems, an underlying asymmetric potential 
acts to break the left-right spatial symmetry. Alterna- 
tive ratchet systems that have recently been investigated 
[5] have symmetric substrates, and rely on ac drives to 
provide symmetry breaking. In each of these cases, the 
particles move through some form of fixed substrate. 

In this work we study a ratchet effect that arises due to 
the collective interactions between driven and non-driven 
particles in a system that does not contain a flxed sub- 
strate. The particular system we consider is a binary 
mixture of interacting repulsive particles which form a 
disordered or glassy media. These particles, referred to 
as the passive species, interact only with other particles, 
and do not respond to an applied ac drive. We introduce 
a small number of a separate species of particles, termed 
the active species, which respond both to the other par- 
ticles and to the ac drive. Using numerical simulations 
we demonstrate that a ratchet effect can be achieved for 
the active particles in this system when the coupling be- 
tween the active and passive particles is strong enough 
and a symmetry breaking is introduced by means of the 
ac drive. This is in contrast to a recent theoretical pro- 
posal of an asymmetric drive ratchet which transports 
undriven superconducting vortices [6]. In addition, for 
specific ac drives, we find current reversals as a function 
of density and temperature. We also show that the mag- 
nitude of the ratchet effect goes through a maximum as 
a function of the number of active particles, due to the 
collective string-like motions of the particles. 

We explicitly demonstrate a ratchet effect for mixtures 
of colloids. Our system is similar to recent studies of 
binary colloidal models driven with ac fields where one 
species of the colloids moves in the direction of the ap- 
plied drive while the other species moves in the oppo- 
site direction [7]. These studies have found interesting 



collective effects in which the moving particles organize 
into lanes. In extensive studies of driven diffusive models 
[8] , where the two species of oppositely moving particles 
are placed on a lattice, a rich variety of noncquilibrium 
phases were observed. Continuum models, with differ- 
ent portions of particles moving in opposite directions, 
have also been studied in the context of pedestrian fiows 
[9] . A remarkable phenomena termed "freezing by heat- 
ing," where the system can jam at higher temperatures, 
has recently been observed in these types of systems [10]. 
Since individual colloids or groups of colloids can be ma- 
nipulated easily with optical traps [11], colloidal systems 
are ideal for experimentally realizing a system in which 
only a fraction of the particles respond to the applied 
drive. The active species could be created using charged 
or magnetic colloids, while the passive particles would 
be charge-neutral or non-magnetic. For example, in re- 
cent experiments, magnetic colloids were driven with ac 
[12] and dc drives [13] through a glassy assembly of other 
non-magnetic colloids. It should also be possible to drive 
a fraction of the colloids using optical or magnetic traps. 

We consider a binary two-dimensional (2D) assembly 
of colloids using the same simulation method as in our 
previous works [14,15]. The overdamped equation of mo- 
tion for colloid i is 



dt 



= tij + (ac + fr , 



(1) 



N 



where ^^j = vtj 
from the other colloids, ^ac 



V^V{r 



is the interaction force 
is the force due to an 
applied driving field, and fr is a thermal force from 
Langevin kicks. The colloids interact via a Yukawa 
or screened Coulomb interaction potential, V{rij) = 
{qiqj/\Ti — rj[) exp(— K[ri — rj|). Here qi is the charge 
on colloid i, 1/k is the screening length, and Tiq-^ is the 
position of particle The system size is measured in 
units of the lattice constant a, which is set to a = for 
most of this work, giving a particle density p = l/a§. We 
take the screening length k = 3/ao. Our system is peri- 
odic in both the x and y directions. The total number of 
particles is = A^^ -f Nb + A^zj, with A'^ colloids of pas- 
sive species A and charge g^, and Nb colloids of passive 
species B and charge qB- We add Nd active colloids 
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FIG. 1. (a) Velocity v vs power-law exponent /3 from Eq. 2 
for /i = 1.0. Ratcheting in the positive direction {v > 0) 
occurs for /3 > 1. Prom top left to bottom left, /2//1 = 1-0 
(flat line), 0.9, 0.8, 0.7, 0.6, and 0.5. (b) v vs (3 for fixed 
/1//2 = 0.5 and, from bottom to top, T2/T1 = 3.25, 2.5, 2, 
1.5, 1.25, and 1. 

with charge qd that respond to an appHed ac drive. For 
the results presented here, we fix qA/qs = 1/2. We have 
also considered other ratios, including the crystalline case 
when qA/qB = ^- To initialize the system, we start at 
a high temperature and anneal cither to T = or to a 
fixed T > 0. The binary particles of species A and B form 
a disordered assembly, with the active colloids scattered 
randomly throughout. After annealing, wo apply a two- 
step ac drive during a single cycle of period r. The first 
stage is a fast drive where a constant force of /i is applied 
for a time ti in the positive x-direction. The second 
stage is a slow drive where a force /2 < /i is applied 
for a time T2 > ti in the negative a;-direction. The ac 
drive has the properties /iri = /2T2, with r = ti +r2, so 
that the particle experiences zero average force during a 
single cycle. The average drift velocity Vd can be obtained 
after many cycles. We conduct a series of simulations for 
different particle densities, system sizes, temperatures, 
qo, ac amplitudes, frequencies, and Nd- 

We first consider a simple model that predicts that a 
ratchet effect can arise in such a system. Recently it was 
shown theoretically and in simulations that, for a single 
colloid driven with a dc drive through a glassy assem- 
bly of other colloids [14], the colloid velocity follows a 
power law ^ oc /" with a = 1.5 for certain parameters. 
The power law appears in regimes where the driven col- 
loid produces collective disturbances of the surrounding 
colloids on length scales of the order of o, the lattice con- 
stant. When the density of the surrounding colloids is 
low, these collective distortions are lost and the veloc- 
ity of the driven colloid is linear with the drive. Linear 
behavior also occurs when the charge qo of the driven 
particle is much smaller than that of the surrounding 
particles. Recent experiments on driven magnetic col- 
loids have also found power law velocity-force relations 
for the driven colloid with exponents a = 1.5 or higher in 
systems with high density [12]. In low density systems, 
the velocity is linear with drive. 

We find that for a group of driven particles in a glassy 
media, the velocity-force curve is also of a power law 
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FIG. 2. Simulation image. Black circles: active particles 
that respond to the ac drive. Open circles: passive particles 
forming the glassy background. Here, Nd/N = 0.016. 

form. In a system with a power law velocity-force rela- 
tionship, ratcheting can occur under the application of 
an asymmetric ac drive. During a single cycle the col- 
loids move a distance Sx = /fn — /I* r2. The net drift 
velocity is then Vd = 5x/{ti + T2)- Using the condition 
/iTi = /2r2 the drift velocity can be expressed as 



Vd 



/2 + /1 



(2) 



If /? = 1, as in the case of a single overdamped particle, 
then = 0, as expected. For any value of /? > 1, how- 
ever, the drift velocity will be positive, > 0. Ratchet- 
ing can still occur if /i is close to /2 as long as /? > 1. 
This model predicts that the ratchet effect becomes more 
efficient for larger values of 0. In Fig. 1(a) we plot Vd vs 
P from Eq. (2), with fi — 1 and ratios of /2//1 = 0.5 
to 1. For /2//1 = 1, = for all (3, as expected. For 
/2 < /i and /? > 1. ratcheting in the positive direction 
(vd > 0) occurs. Vd increases with larger [3 and larger 
f'l/fi- The velocity reverses for /3 < 1. This is an un- 
physical limit for the colloid system, where the lowest 
possible value is P — I. However, some physical systems 
such as non-Newtonian fluids could have (3 < 1. 

Since the velocity becomes linear with the drive at low 
densities or high drives, it should be possible to select 
ac drive parameters that reverse the ratchet effect and 
give current reversals for f3 > I. This can be achieved by 
setting n < T2, where the inequality is not too large. In 
the highly nonlinear regime the particle will still ratchet 
in the positive direction; however, upon approaching the 
linear regime, the particle will ratchet in the negative 
direction. In Fig. 1(b) we show v vs /3 from Eqn. (2) for 
/2//1 = 0.5 and increasing t^/ti, where a crossover from 
a positive to negative ratchet effect occurs at /3c. Here 
Pc increases with increasing tiIt\. 

We next consider simulation results. In Fig. 2 we show 
a real space image of the simulated system with Nu = 11 
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FIG. 3. Particle velocity v vs time t/r. Dotted line: v for 
an isolated overdamped particle under an asymmetric ac drive 
with /i = 0.56, /2//1 = 0.25, and ri/r2 = 0.25. Solid line: v 
in a system of interacting particles. Inset: Drift velocity va vs 
t/r for (dotted line) a single overdamped particle, and (solid 
line) the interacting particle system. 

active particles (black circles) and Na + Nb — 678 pas- 
sive particles (open circles) after many cycles of the ac 
drive. The active particles, which were originally scat- 
tered throughout the system, have clustered into a lane, 
and they also drift over time in the positive x-direction. 

In Fig. 3 we plot the velocity vs time of a single over- 
damped colloid in the absence of any other particles 
(dashed curve), showing short large positive pulses and 
long small negative pulses of velocity. Here we consider 
the case of /2//1 = 0.25 with /i = 0.56. The average ve- 
locity of the isolated colloid is zero in a single period. The 
solid line shows the velocity of a driven colloid in the in- 
teracting particle system with = 370 and No = 11, in 
a regime where the dc velocity-force curves have a power- 
law form. Here, the magnitude of the particle velocity v 
during both the positive and the negative portions of the 
ac drive cycle is smaller than that of the single particle, 
and shows considerable fluctuations. The integrated ve- 
locity during the positive portion of the ac cycle is on 
average larger in magnitude than the integrated veloc- 
ity during the negative portion of the cycle. To show 
this more clearly, in the inset of Fig. 3 we plot the time- 
averaged drift velocity for each system over many ac 
cycles. For the single particle case (lower dashed curve), 
frf = 0, while for the interacting particle system, Vd > 0, 
with an long time average drift velocity of < Vd >~ 0.03, 
indicating that a positive ratchet effect is occurring. 

We next consider the dependence of the ratchet effect 
on various system parameters. In Fig. 4(a) we plot the 
drift velocity Vd vs the fraction of driven particles, Njj /N, 
for fixed ac drive parameters /i = 0.56, /2//1 ~ 0.25, and 
ti/t2 = 0.25, with a/ao — 1. Vd goes through a 




FIG. 4. (a-c) Drift velocity Vd obtained from smiulations 
with /iTi = f2T2 and /2//1 = 0.25. (a) Vd vs the fraction 
of driven particles Nd/N for a density of p — l/af). (b) 
Vd vs particle density p for fixed Nd /N = 0.03. (c) Vd vs 
temperature T/T^ for p ^ l/ag and Nd/N = 0.03. (d-e) 
Ratchet reversals for a system with fiTi/f2T2 = 0.9: (d) va 
vs density p. (e) Vd vs temperature T/Tm for p = 1.26/ao. 

maximum near No /N = 0.03, and then slowly falls off for 
higher Njj/N before reaching Vd — just below N^/N = 
1.0. We have also considered different system sizes by 
varying a/aQ. We obtain the same behaviors, indicating 
that the ratcheting is not a finite size effect. It is simple to 
understand why Vd ui the limit of N^/N — > 1, since 
here all the particles move in unison with the ac drive, 
and the plasticity responsible for the nonlinear velocity- 
force response is lost. For a single driven particle No — 
1, we found a dc velocity- force scaling exponent oi [3 — 
1.5. As the number of driven particles No increases, (3 
increases to (3 — 2 near No/N — 0.03, and then decreases 
back toward /3 — 1 as No/N increases further. This 
result is consistent with the predictions of Eq. 2, where 
a larger ratcheting effect is expected for larger exponents 
/3. 

In Fig. 4(b) we consider a system with fixed No/N = 
0.03 and ac drive parameters, and plot Vd vs the system 
density p measured in units of I/oq. For low densities 
p < 0.5, the active particles interact only weakly with 
the passive particle background, and the dc velocity-force 
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curves are close to linear (/? ?a 1) so little ratcheting 
occurs. For large densities p > 1.3, the ratcheting drift 
velocity Vd decreases since it becomes more difficult for 
the active particles to pass the passive particles. Wc 
have also fixed the system density and varied the charge 
qd of the active particles. The behavior is similar to 
Fig. 4(b): the ratchet effect is lost at small qD when 
the velocity-force curve becomes linear. In Fig. 4(c) we 
show the effect of a nonzero temperature on the ratchet 
effect. We plot Vd vs T/Tm, where Tm is defined as the 
temperature at which a completely ordered single-species 
colloid assembly with charge qn and density p = I/oq 
melts. For high T, the velocity increases linearly with 
force, so the ratchet effect is lost. 

Equation 2 predicts that a reverse ratchet effect should 
occur for /? < 1, as shown in Fig. 1(a), or for /iTi ^ f2T2- 
The reverse ratchet effect can occur in systems where 
the media responds easily to slow moving particles but 
becomes stiffer under faster perturbations. For example, 
in granular media the grains respond in a fluid manner 
to a slow moving object, while for a fast moving object 
the granular material jams and acts like a solid. Many 
polymer systems show a similar response. In order for a 
ratchet effect to occur in a substrate-free system, there 
must be a combination of the nonlinear behavior of the 
particle motion and the left-right symmetry breaking by 
the ac drive. In general the ratchet effect we observe 
should occur for any media that has a nonlinear viscosity- 
frequency response. 

In Fig. 4(d,e) we consider two cases where flux rever- 
sals occur with /1T1//2T2 = 0.9. For the parameters cho- 
sen here, the drift velocity of an isolated particle would 
be Vd = —0.016. In Fig. 4(d) we plot Vd vs density p for 
Nd/N — 0.03. For low density p < 0.3, the velocity- force 
response curve is linear (/3 = 1) and the particle moves in 
the negative direction as expected for this ac drive. As 
the density increases, the velocity- force response curve 
becomes nonlinear, and Vd crosses over to a positive value 

> as the ratcheting effect reaches a maximum near 
p = l/ftQ. The drift velocity drops at high density just 
as in Fig. 4(b). In Fig. 4(e) we show V vs T/T^ for the 
same system with a fixed density oi p = 1.26/aQ. For low 
temperatures, the velocity-force response curve is non- 
linear, and a ratchet effect with positive velocity occurs. 
As the temperature increases, the drift velocity Vd drops 
and crosses back to a negative value as the system enters 
the linear response regime. Vd saturates to Vd ~ —0.01, 
slightly smaller in magnitude than the value expected 
for a single isolated particle. We have also considered 
the effect of changing at fixed density and T = 0, and 
observe behavior similar to that shown in Fig. 4(d). 

In conclusion we have investigated a "glassy ratchet" 
effect, which occurs in a glassy system of interacting par- 
ticles without a quenched substrate. We considered the 
specific case of a disordered assembly of colloids with 
overdamped dynamics where only one colloid species cou- 



ples to an ac drive. The ac drive is asymmetric, with 
a large positive drive applied for a short time followed 
by a small negative drive applied for a longer time, and 
obeys the constraint that the net applied force on the 
driven colloids is zero in a single cycle. A ratchet effect 
occurs when the driven and the non-driven colloids are 
coupled, producing a nonlinear velocity-force response. 
We find that the dc drift velocity has a maximum value 
as a function of the fraction of driven particles, tempera- 
ture, and particle density. A simple model for the ratchet 
effect agrees well with our simulation results, including 
the occurrence of ratchet reversal regimes. The system 
considered here can be realized experimentally for mag- 
netic or charged colloids driven through non-magnetic or 
uncharged colloids. Our results may also have applica- 
tions for new types of electrophoresis devices. 
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